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Abstract. We present recent results on the ^-dependence of four-dimensional SU(N) gauge 
theories, where 6 is the coefBcient of the CP-violating topological term in the Lagrangian. 
In particular, we study the scaling behavior of these theories, by Monte Carlo simulations at 
imaginary 9. The numerical results provide good evidence of scaling in the continuum limit. 
The imaginary 6 dependence of the ground-state energy turns out to be well described by the 
first few terms of related expansions around ^ = 0, providing accurate estimates of the first few 
coefficients, up to 0(6''). 



Four-dimensional SU{N) gauge theories have a nontrivial dependence on the parameter 6 
which appears in the Euchdean Lagrangian as 

Ce = {l/4)F^,{x)F^,ix)-ieq{x), q{x) = /{QATT^)e^,p^F%{x)F''p„{x), (1) 

where q{x) is the topological charge density. The ground-state energy density F{9) behaves as 

T{e) ^ F{e) - F(0) = (1/2) xO^s{0), (2) 

where x is the topological susceptibility at = 0, 

X = Jd^x{q{x)q{0))e=o = {Q^)e=o/V, Q = Jd'^x q{x), (3) 

V is the spacetime volume and s{0) is a dimensionless even function of such that s(0) = 1. 
Assuming analyticity at ^ = 0, s{0) can be expanded as: s{9) = 1 + 62^^ + ^4^^ + ■ ■ ■ > where 
only even powers of 6 appear. Large- scaling arguments applied to the Lagrangian (1) indicate 
that the relevant scaling variable in the large- limit is ^ = 0/N. This implies that in this limit 
X = 0(1), while the coefficients 621 are suppressed by powers of N, i.e. 621 = 0(A"~^*). 

Due to the nonperturbative nature of the 9 dependence, quantitative assessments have largely 
focused on the lattice formulation of the theory, using Monte Carlo (MC) simulations. However, 
the complex nature of the term in the Euclidean Lagrangian prohibits a direct MC simulation at 
0. Information on the 9 dependence of physically relevant quantities, such as the ground state 
energy and the spectrum, has been obtained by computing the coefficients of the corresponding 
expansions around 9 = 0. The coefficients of s{9) can be determined from appropriate zero- 
momentum correlation functions of q{x) at 9 = 0, which are related to the moments of the ^ = 
probability distribution P{Q) of the topological charge Q. Indeed 
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Figure 1. Distribution of the ratio 
Qt/tti for /3 = 6.2 configurations at 
9l = and Ol = 8.7 (Oi ^ 1.5). 



etc. They parameterize the deviations of P{Q) from a simple Gaussian behavior. It has 
been shown that correlations involving multiple insertions of the topological charge can 
be defined in a nonambiguous, regularization independent way, and therefore 621 are well 
defined renormalization group invariant quantities. The numerical evidence for a nontrivial 
6'-dependence, obtained through MC simulations, appears quite robust. We refer the reader to 
Ref. [1] for a recent review. On the other hand, MC simulations at = have only made it 
possible to estimate the ground-state energy up to 0{9^). The large- prediction 62 = 0{N^'^) 
has been already supported by numerical results; the calculation of the higher-order terms would 
provide a further check of the power-law suppression predicted by large- arguments. 

In this paper we consider imaginary values of 9, which make the Euclidean Lagrangian (1) 
real, thus making MC simulations possible. For further details and references, the reader may 
consult our Ref. [2]. Assuming analyticity at ^ = 0, the results provide quantitative information 
on the expansion around 9 = 0. Indeed, fits of the data to polynomials of imaginary 9 may 
provide more accurate estimates of the coefficients, overcoming the rapid increase of statistical 
errors observed at 9 = 0. Perturbative renormalization-group (RG) arguments indicate that 9 is 
a RG invariant parameter of the theory, thus the continuum limit should be approached while 
keeping 9 fixed to any complex value. We find that this is indeed supported by the numerical 
data for the 4D SU(3) lattice gauge theory, at least for \9\ < vr, which are well described by the 
first few nontrivial terms of the expansion around 9 = 0. 

Introducing the real parameter 9i, defined hy 9 = —i9i, Eq. (2) leads to: 
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The nonperturbative formulation of the above theory on the lattice requires a discretization 
of the action, Sl — 9lQl', for Sl we use the plaquette gluon action, while for Ql we employ 
the "twisted double plaquette" operator {Ql = J2x QLix)). Notice that this is not the only 
possible choice for qi ; the only requirement is that it have the correct continuum limit when 
a — )• (a: lattice spacing). In the continuum limit qiix), being a local operator, behaves as 

qL{x) ^a^Z,q{x) + 0{a^), (7) 
where Zq is a finite function of the bare coupling qq, going to one in the limit /3 = 2N/gQ 



00. Thus, we have the correspondence: 9i 



Zq9L, 



apart from 0{a ) corrections. The 



renormalization Zq may be evaluated by MC simulation at ^ = 0, computing 

Zq = {QQL)e=o/ (Q^)e=o , 



(8) 



where Q is an estimator such as those obtained by the overlap method or the cooling method, 
which are not affected by renormalizations, nor by nonphysical contact terms. Thus, the ratios 



{Q)eJ{Q^)e=o = 9^- 2b29f + %^9f + {Q%/{Q^)e=o = 1 - 662^^ + 1564^.^ + (9) 
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Figure 2. vs 0, = Z^Ol- Figure 3. {Q)eJ{Q%L vs 0^ = Z^Ol- 

are expected to have a well defined continuum limit as functions of 9i . 

We have carried out MC simulations of the 4D SU(3) lattice gauge theory, at /3 = 5.9, 6, 6.2, 
for lattice sizes L = 16, 16, 20, respectively; the simulations are carried out both at 6*^ = 
and 0L 7^ 0, within the region \6i\ < vr. Since our numerical study requires high-statistics 
MC simulations, we choose the cooling method as estimator of the topological charge Q. The 
topological charge has been measured on cooled configurations (by locally minimizing the lattice 
action), using the twisted double plaquette operator. As is well known, this procedure leads to 
values Qt — ka, where k is an integer and a ^ 1. Once we determine a, we assign to Q the 
integer closest to Qt/ct. This cooling method for estimating Q, though less rigorous than the 
significantly more expensive overlap method, produces results in good agreement with it. 

Fig. 1 shows the distributions of the ratio Qt/a of /3 = 6.2 cooled configurations at = 
and 6l = 8.7 {9i = ZgOi 1.5). We note that these distributions cluster around integer values, 
also for rather large values of Qt/a, both for 9^ = and 9^ = 8.7. 

MC simulations at 9 = were performed at /3 = 5.9, 6.0, 6.2. Over 40 million sweeps per 
value of /3 were produced. The results for X) ^2, ^4 and Zg are reported in Ref. [2]. Providing 
our improved estimate of high-order coefficients, such as 64, turned out to be very hard in 
9 = MC simulations, requiring huge statistics. The results for 64 are consistent with zero, 
suggesting the bound [64 1 < 0.005, which is improved in 7^ runs. Our estimates of Z„ 
Zq{l3 = 6.0) = 0.135(1) ) reduce the uncertainty on Zg as produced by other methods. 

MC simulations at 9^ ^ are slower by approximately a factor of three, due to the complexity 
of the action. In 9l ^ runs, ^ 3 million sweeps were produced for each value of f3 and 9l. 
Figs. 2 and 3 show results for the ratios ((5)6»i/(Q^)e=o and (Q)6»i,/(Q^)0^ ) versus 9i = Zq9L (cf. 
Eq. (9) ). The MC data at different /3 values follow the same curve, providing evidence of scaling. 
Scaling corrections, expected to be 0{a?), are quite small, and tend to increase with increasing 
9i. This good scaling behavior corroborates the existence of a nontrivial continuum limit for 
any value of 9i. Fitting our data to Eqs. (5,6,9) improves significantly the 6* = results. In 
particular, a much smaller bound on 64 is obtained: [64 1 < 0.001 ; also, we find: 62 = —0.026(3), 
which is clearly more precise than the estimate obtained from 9 = runs only: 62 = —0.029(7). 

Besides allowing more precise determinations of the 9 expansion coefficients of the ground- 
state energy and other observables, using imaginary 9 values might turn out useful in overcoming 
the dramatic critical slowing down of topological modes, by performing parallel tempering 
simulations with a set of 9 values including 0=0; this is an exact MC algorithm for the model. 
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